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The stored ultra-cold neutrons have been developed. A high density ultra-cold neutron 
gas has been recently produced by using the nuclear spallation method. We investigate the 
thermodynamic properties of the quantum ultra-cold neutron gas in the Earth's gravitational 
field. We find that the quantum effects increase temperature dependence of the chemical 
potential and the internal energy in the low temperature region. The density distribution of 
quantum ultra-cold neutron gas is modified by the Earth's gravitational field. 

§1. Introduction 



On 

■ Neutrons are difficult to trap in any container because of their neutral nature. 

Very slow neutrons, however, are found to be reflected by a metallic wall, and there- 
fore, can be stored in metallic containers; see Ref. [Tj) and references therein. They 
are the so-called Ultra-cold Neutrons (UCNs). The technology of producing UCN 
has been developed in order to search for the time reversal violating electric dipole 
moment of the neutron and for the most accurate measurement of the neutron j3 
decay lifetime. The UCN scattering is applied to the study of condensed matter 
properties [2]). Recent experiments have succeeded in creating UCN gas by using 
nuclear spallation method, where the neutron density is not limited by Liouville's 
theorem ED ,|H), ED, ED. The density of UCNs coming from a nuclear reactor never in- 
creases under this restriction. Therefore the nuclear spallation method is expected 
' to create high density UCN gas and also to study quantum neutron systems. 

Turning to nuclear physics, neutron systems have been studied to discover inter- 
esting phenomena originating in the strong interaction like the neutron starEJ. This 

. (=10 39 cm -3 ). In the current experimental state ED the UCN gas is very dilute, with 

density ~ 10 2 cm -3 , At such low densities, thermodynamic properties are described 
by the Maxwell-Boltzmann distribution: the mean separation length between neu- 
trons is R = (density) -1 / 3 ~ 2 x 10 -1 cm, the thermal de Broglie wave-length is 
A = hf '(SmfcfiT) 1 / 2 ~ 8 x 10 -6 cm at 1 mK, where m is neutron mass, ks is the 
Boltzmann constant and T is the temperature. Due to the very long wave-length, 
the UCN can be reflected and stored in metallic containers. 

In this paper, we investigate the thermodynamic properties of the quantum 
neutron gas in the Earth's gravitational field. With developing UCN experiments 
to create high density or low temperature gas, a treatment of quantum statistical 
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mechanics is needed. The neutron density is expected to increase under the influence 
of the gravity, which is even more enhanced as the temperature becomes lower. There 
the chemical potential and the internal energy increase as compared to the case of 
the free Fermi gas. The density is the highest at the bottom of the container, and it 
becomes smaller with increasing height and with increasing temperature. In Sec. [2] 
we show the wave function and the eigenenergies of a neutron trapped by the Earth's 
gravitational field. Sees. [3] and [3] are devoted to the internal energy and the density 
distribution by using the chemical potential. Finally, we conclude the present study 
in Sec. 



§2. Neutron in the Earth's Gravitational Field 



Let us first consider the one-neutron wave function trapped by the Earth's grav- 
itational field. The Shrodinger equation reads 

h 2 d 2 it>(z) , „ , , , . 

where g is the standard gravitational acceleration, the z-axis is along with the 
Earth's gravitational field and z = is the bottom of the container. By introducing 
a = 2m 2 g/h 2 , e = (a/mg)E z and £ = (az — e)/a 2//3 , Eq. (|2-1|) is rewritten as a 
dimensionless Schrodinger equation 

^-«(0-O. (2.2) 

The general solutions are given by the well-known Airy functions (see Ref. [8])) 

<K0 = N a Ai(0 + N b Bi(£) , (2-3) 

where N a and iVj determine the weight of the two independent functions. Assuming 
the boundary condition that lim^oo ip(z) = 0, the factor of Nb should be zero, 
because Bi(£) diverges at infinity. Therefore the wave function is proportional to 

^ z )= Na Ai[^^ . (2-4) 
By imposing another boundary condition ip nz (0) = 0, the wave function must obey 



Ai {-Jf-3)= > ^ 

which determines the eigenenergies because the neutron is reflected by the bottom 
of the container. By using the asymptotic behavior of zeros of the Airy functions [8]) , 
the eigenenergies q| are given by 

En, - (^f) V3 (f ) (4n, -!)■/.,„,_!, 2 ... . (2-6) 



This is a good approximation at 1 % or better 
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Substituting g = into Eq. ()2-6p . obviously, the energies are vanished. This reflects 
the fact that the z-direction is not bounded. The normalization N a in Eq. (|2-4p is 
determined from the orthonormalization condition of the wave functions, 

1pn(z)ll>m( z )dz = $n,m ■ (2-7) 

By using the integration formula from Ref. ED, the normalized wave functions are 
given by 

*- (z) = m^j^)\ M rs^j ■ (2 ' 8) 

Furthermore the very slow neutron is reflected by a metallic wall in x- and y-axes. 
For simplicity we consider a square- well potential of infinite height. Thus the energies 
are E n = 7t 2 h 2 n 2 /2mL 2 , where L is the size of square- well and n is a positive integer. 
Therefore, the total energies are written by the sum with respect to the three-axes: 

where n x , n y and n z are arbitrary positive integers. 

§3. Density of States, Chemical Potential and Internal Energy 

To derive the thermodynamic properties, let us calculate the density of states 
by using Eq. {232): 



dn x dn y dn z 2m ( L \ 2 1 / 2 / E Uz 



p(E nx ,E ny ,E nz ) - ^ ^ ^ - h2 j ^ E ^ . (3-1) 

The chemical potential /i can be determined from the following condition: 

f°° p(E n ,E n ,E n ) 

N = 2 l ^.^^. eW £,;x,X-^ + i - (3 ' 2) 

where N is the total number of particles and the factor of 2 comes from spin degen- 
eracy. The coefficient /3 is defined as (A^T) -1 , where ks is the Boltzmann constant 
and T is the temperature. In the limit T — > the Fermi energy is 

h 2 (l^ 2 m 2 gN\ 2/b 



€F 2m 



W^T 2 ) • 



Notice that it is different from that of the free Fermi gas; the volume depends on the 
Earth's gravitational field; the power in Eq. (|3-3p is 2/5, which is related with the 
density. To simplify the density of states, putting E Ur = E Ux + E ny (n x = n r cos 9 
and n y = n r s'm9) in Eq. (|3-2p . it is rewritten down 



8e^ 
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where we use e F and 



o J o 



dE n JE ny 2h 2 ir 2 f°° f°° r/\ JQ [°° 2h 2 7r 2 n r [°° 

" " dn x dn y = / do j ^ — an r = ir / dh nr 



\J E nx E Tly mL 2 J J J J mL 2 

(3-5) 

The density of states can be redefined as 

15JV 



p(E nr ,E n ,) = —^E^. (3-6) 
8e F 

Setting 7] = /3fi, £ = j3E nz and v = /3E nr for Eq. (|3-4p . f3 is regarded as a function of 
V, 



15 r r c i/2 ^ V /5 ( 5 r c 3/2 ^V /5 

t/ io if^F+i^H =(2/0 ^n dc ) 



where the integration over v is performed by integrated by parts. Therefore the 
chemical potential is given by 



H 7] 



-2/5 

nil I tt^-t^C I , (3-8) 



5 roo ^3/2 



e F f3e F '\2J e^ + l 

which is normalized by e F . In the low temperature approximation, the integration 
of Eq. (GEED is analyzed in Refs. fTU]) .[TT )) . 

jH ^Sri^ = l^ /2 + T^ 1/2 + • • • > for e-" « 1 . (3-9) 
Thus the chemical potential can be calculated as 

ii = l-^ + ,fore-^«l. 

£F 2 V e F J 

The internal energy is given by 

= w f 1 r + r r f 2v , «* 



(3-10) 



(3-11) 



In Fig. (H w e plot the temperature dependence of the chemical potential and the 
internal energy using Eqs. (j3-7U3-8l 13-lip . which are compared with those of the free 
Fermi gas. The temperature dependence of the chemical potential is 6 times larger 
than that of the free Fermi gas with the same Fermi energy in the low temperature 
region. The internal energy also has strong temperature dependence at the low tem- 
perature region. As in the case of the classical gas like atmosphere, this means that 
the specific heat of quantum UCN gas increases by the Earth's gravitational field. 
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Fig. 1. The chemical potential and the internal energy as functions of the temperatures: the solid 
lines present the results of Fermi gas under the Earth's gravitational field, the dashed lines are 
those of the free Fermi gas. 



§4. Density Distribution 

It is difficult to observe, for example, the specific heat of UCN gas because the ex- 
ternal parameters of pressure or temperature are not easily changed in experiments. 
However the density distribution can provide information on the thermodynamic 
properties, and it can be observed by neutron counting experiments. To clearly show 
the z- and T- dependencies of the density distribution, following an example shown 
in Ref I12p . the condition of chemical potential in Eq. (|3-2p can be rewritten as 

N _ 1 /"^3 1 _ f°° P(En r ,E n J j 

~ (2TThf J P eif/lm+rngz-ri/kBT + \~ J Q e /B(E nr +E nz -„) + ^n^n. ■ 

(4-1) 

Here we define the single particle density as 

1 r 3 i 

n(T ' Z) = (2^3 J d P eiiP/lm+mgz-ri/kBT + { ' ^ 



then Eq. (|4-ip becomes 



N= I dx dy I dzn(T,z). (4-3) 



o 



Eq. ()4-2p can be calculated in the limit T — > as 

(2m(ep — mqz)) 3 / 2 „ A . 
n (°' Z ) = ^ Q^J > fOT e F > m 9 Z > ( 44 ) 

while it becomes zero for tp < mgz. Similarly, putting x = f3p 2 /2m, the T-dependence 
of the density distribution at z = is given by 

oo 1/2 



n(T,0) _ 3, _ m _ 3/2 f°° X 
n(0,0) 



3 f°° x 1/z 

-(e F /3y 3/2 / — * d X , (4-5) 
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which is normalized by n(0, 0). In the low temperature approximation, the integration 
is analyzed in Refs. HQ ]) . ITT]) . 

/ d X = -v V2 + — V 1/2 ■ ■ ■ , for e-" « 1 . (4-6) 

7 eX-v + 1 A 3 ' 12 ' ' v ; 

Thus the density distribution at z = is calculated as 

n(r ' 0) -l-^f^V + ...,for e -/^«l, (4-7) 



n(0,0) 8 V e F J 

by using the definition of f3 in Eqs. (|3-7| 13 -9|) . In Fig. [21 we plot the density dis- 
tribution n(T,z) using Eqs. (|3-7l !4-2p . The density distribution with respect to z 
coordinate rapidly drops with mgzjep in the low temperature region, and it also 
becomes flat in the high temperature region. Finally, we discuss the properties of 




Fig. 2. The density distribution as functions of the temperature ksT/eF and mgz/tF 

quantum UCN gas. In Fig. [3l we plot the density at the bottom of the container 
as a function of epfk^ in terms of the double logarithmic form. For example, when 
the Fermi energy of quantum UCN gas is €fI^b ~ 1 mK, the density becomes 
~ 10 16 cnr 3 at the bottom and the height also becomes ep/mg ~ 8 x 10 1 cm. In 
fact, the mean separation length R ~ 5 x 10 -6 cm is close to the thermal de Broglie 
wave-length A ~ 8 x 10 -6 cm. 



§5. Conclusions 



The previous sections provide general thermodynamics properties and the den- 
sity distribution of the neutron gas under the Earth's gravitational field. The chemi- 
cal potential and the internal energy have strong temperature dependence at the low 
temperatures because the density of sates depends on the square root of the eigen- 
energies in the Earth' gravitational field. Thus it can be expected that the variation 
of specific heat with temperature is larger than that of free Fermi gas. Instead of the 
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l°9lO E F /k B ( K ) 

Fig. 3. The density at the bottom of the container as a function of the Fermi energy cf/^b 

thermodynamic indexes, the density distribution shows how the quantum UCN gas 
is affected by the Earth's gravitational field. The density is the highest at the bottom 
of the container and it becomes small with increasing mgz/eF and temperature T. 
Assuming that the average kinetic energy of UCNs becomes much smaller than the 
Fermi energy, the quantum UCN gas would be created. 
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